COBE Observations of Interstellar Line Emission 
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Abstract. The FIRAS instrument on the COBE satelhte has conducted 
an unbiased survey of the far-infrared emission from our Galaxy. The first 
Q ■ results of this survey were reported by Wright et al. (1991), and later 

O , conclusions were reported in Bennett et al. (1994). I report the results of 

Q ' analyses of this spectral survey, which includes emission lines from 158 ^um 

C+, 122 /xm and 205 /im N+, 370 ^m and 609 /xm C°, and CO J=2-l 
through 5-4. The morphological distribution along the galactic plane 
(6 = 0°) for all the spectral line emission, and the high galactic latitude 
intensities of the strong and 205 /im N+ emission are discussed. The 
high galactic latitude intensity of the 158 /xm fine structure transition 
^ ; from C+ is (C+ 158 /xm) ^ (1.43 ± 0.12) x 10"^ esc \b\ erg cm'^ s'^ sr"! 

for \b\ > 15°, and it decreases more rapidly than the far infrared intensity 
CN I with increasing galactic latitude. 



1. Introduction 



O I The interstellar medium is heated primarily by photons interacting with grains, 

-4— > ' producing hot electrons by the photoelectric effect (Watson 1972). Other heating 

I is supplied by the excess energy given to electrons detached from atoms during 

photoionization, and by the excess energy given to atoms during the photodis- 
sociation of molecules. The hot electrons and atoms collide with the thermal 
^ ' particles and add energy. To maintain a steady state, this added energy must 



be radiated by the gas. The temperature of the gas will increase until the radi- 
ation process with lowest excitation threshold becomes effective. In the neutral 
atomic interstellar medium, the fine structure line of the C"*" ground state with 
an excitation energy corresponding to 91 K is the first effective cooling chan- 
nel to open, so the C"*" line at 157.7 /xm is expected to be strong (Dalgarno & 
McCray 1972). The ionization threshold for C is 11.26 eV, so it will be ionized 
even when H is neutral. O and N have higher ionization thresholds and will 
be neutral in H I regions. The ground electronic state of O is an inverted fine 
structure triplet, so the first transition from the ground state corresponds to the 
63 fxm line, and thus requires a considerably higher temperature for excitation. 

In the warm ionized medium, the CNO elements will all be singly ionized. 
The cooling channel with the lowest excitation energy will be the N"*" line at 
205 /xm, but the 158 ^m line of C"*" will still be strong. With a typical tempera- 
ture of a few thousand degrees in the warm ionized medium, there will be many 
more channels for cooling, so none of these lines will dominate the cooling in the 
way the C"*" line dominates in H I regions. 
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Figure 1. FIRAS Right High Spectra of the Galactic Center and the 
Lockman Hole. 

In molecular gas, the rotational transitions of the CO molecule provide low 
threshold energy cooling channels. The density and radiation field in regions 
with CO favor the presence of neutral C instead of C"*". The ground electronic 
state C*^ is a fine structure triplet with lines at 609 and 370 fiia. 

In this paper I will discuss the observations of these lines by the FIRAS 
instrument on the COBE satellite. These data were first discussed by Wright 
et al. (1991) which included the first detection ever of the 205 /um N"*" line. 
Bennett et al. (1994) discuss the distribution of these lines around the sky. 

2. Observations 

The FIRAS instrument is a polarizing Michelson interferometer first described 
by Martin & Puplett (1970). Using free-standing wire grids as polarizers allows 
a single instrument to measure many octaves of spectrum at once. Figure |l| 
shows spectra of the Galactic Center pixel and the pixel with the lowest Nh 
(Lockman, Jahoda and McCammon 1986) taken from the Right High channel 
of FIRAS. The high frequency part of the spectrum (> 20 cm~^) was separated 
from the cosmic signal by a dichroic beamsplitter, but the residual reflectance 
of this filter allows one to see most of cosmic blackbody spectrum as well. 

The FIRAS instrument uses bolometric detectors, and is thus approximately 
equally sensitive to lines at any frequency. However, variations of filter trans- 
mission and modulation efficiency modify this equality and must be determined 
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using calibration observations. The in-flight calibrators are essentially black- 
bodies, while the galactic spectrum is a very dilute modified Planck function. 
Thus the 50 part per million precision achieved in the search for deviations of 
the CMB from a blackbody (Fixsen et al. 1996) does not imply 50 part per 
million precision in the galactic dust and line fluxes. In particular, small ripples 
in the gain calibration could be confused with weak line fluxes in regions where 
the continuum is high. 

The bright and N+ lines are easily visible in the spectrum of the galactic 
center, and by fitting a line profile plus baseline to the spectrum of each pixel, 
maps can be generated. Figure |2| shows a map of the C"*" line flux. 

Higher SNR on the galactic signal can be obtained by averaging all of the 
pixels together. But an arbitrary isotropic signal must be included in the fit, as 
well as the dipole variation of the CMB temperature. The following approach 
was used: 

1. Make maps of the cosmic temperature variation and the galactic dust 
distribution using an initial guess for the mean galactic spectrum, Go{v)'- 

Ul, b) ^ B,(To + AT(/, b)) + g{l, b)Go{u) 

2. Do a spectral fit to derive the isotropic cosmic spectrum, Io{v)] the spec- 
trum of the dipole, D{i'); and the mean spectrum of the galaxy, G{v): 

b) « Io{y) + D{v) cos e + G{v)g{l, b) 

where is the angle between the LOS (/, b) and the hot pole of the dipole. 
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Figure 3. FIRAS dust map 6). 

3. Iterate using G(z^) instead of Go(i^). 

This procedure can be applied to a subset of the sky as well. Note that changes 
in the isotropic spectrum have no effect of the derived G{u), which is determined 
solely from the variations of the spectrum with position on the sky. 

Wright et al. (1991) applied this procedure to the early data from FIRAS, 
and obtained a map of g{l,b) over part of the sky and a spectrum G{i') which 
was analyzed for line strengths and dust emission parameters. Figure |^ shows 
g{l, b) for the entire FIRAS data set. Note that the polar flux (slope with respect 
to esc is 0.05 units on this map. 

The spectrum G{iy) derived only from data taken away from the galactic 
plane is given by Fixsen et al. (1996). Figure Q shows this G{u) after a modified 
Planck function is subtracted. The 3-2 and 4-3 CO lines have about the same 
strength in this spectrum as they do in the all-sky G{u), while the 5-4 line is 
stronger. However, the SNR of this spectrum is quite low since all of the strong 
galactic signal has been ignored. 

3. CO Analysis 

FIRAS detected the CO lines from 2-1 to 5-4. In terms of power in erg/sec, the 
5-4 line is quite strong. The excitation of CO to the J = 5 level which has a 
radiative lifetime of only one day would seem to be quite difficult in the high 
galactic latitude molecular clouds seen by Magnani, Blitz and Mundy (1985). 
This excitation looks more reasonable if the CO line strengths are plotted in 
units of K-km/sec, as shown as in Figure ^. The almost perfect fit of these 
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Figure 4. Fixsen et al. (1996) G{v) minus a modified blackbody. 
The positions of the 2-1 through 5-4 CO hues are marked. 

line strengths to exp(— 6Jf/) (dashed line) suggests a thermal distribution, but 
since the energy of the CO levels is cx J(J -|- 1) a mixture of temperatures is 
needed: the solid curve shows a model of optically thin CO with 88% at 5 K 
and 12% at 19 K. While this gives a good fit, the number of parameters equals 
the number of detected lines, so there are no remaining degrees of freedom. 
This two temperature model predicts a CO 1-0 line flux of 0.15 K-km/sec at 
the galactic pole. The standard ratio of CO 1-0 line flux to molecular hydrogen 
column density, N{H2)/I{CO) = 3 x W^/^cth^ K km/sec) (Young & Scoville 
1991), gives N{H2) = 5 x 10"^^ so this molecular gas is about 1/3 of the atomic 
gas at high latitudes. Blitz, Bazell &: Desert (1990) estimate that the mass of 
the high latitude molecular clouds is 33% of the neutral atomic mass, implying 
that the FIRAS CO flux is due to the smeared out emission from the MEM 
clouds. The high latitude CO clouds have lower densities and lower extinctions 
than most molecular clouds, so a fluorescent cascade following vibrational or 
electronic excitation could provide most of the high- J CO. 

4. Conclusion 

The FIRAS data on the Milky Way show the primary cooling mechanisms in the 
ISM. 99.7% of the power is radiated by dust heated directly by the interstellar 
radiation field. Gas in the ISM radiates 0.3% of the power in the C"*" line at 
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Figure 5. FIRAS CO line strengths at g{l,b) = 1. The sohd hne 
shows a two-temperature model fit. 

158 /xm. The next most luminous species seen by FIRAS is N"'", with 0.03% of 
the power in the 205 jim line. C° and CO each radiate about 0.003% of the total 
power. 

Acknowledgments. This work was made possible by the dedication of 
hundreds of scientists and engineers who worked to make the COBE project a 
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Discussion 

Dr. Nakagawa: What is the temperature range derived from the CO lines? 
Prof. Wright: The fit gives 5 K and 19 K. 

Prof. Greenberg: When you say 88% of the CO at 5 K and 12% at 19 K do 
you imply that the CO is thermally coupled to the dust? I find this difficult to 
believe. Is this just a coincidence? 

Prof. Wright: I hope so! Collisional heat transfer between gas and dust is quite 
slow. 

Dr. Taniguchi: Your continuTim and [C II] maps show an emission component 
North of the Galactic Center. What is it? 

Prof. Wright: That is the Sco-Oph dark cloud complex. 
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Dr. Taniguchi: Is the northern extension at / = and the southern extension 
at the anti-center due to warping of the disk? 

Prof. Wright: No, those are both very nearby cloud complexes. The warp of the 
disk is a much smaller angle. 

Dr. Giard: What is the upper limit on the ortho:para ratio in H2O based on 
the apparent absorption at the Galactic Center? 

Prof. Wright: First, we do not claim the absorption at / = is real. But if it is, 
the line strengths shown in Figure 2 of Bennett et al. (1994) imply an equivalent 
width of Wp = (3 lb 4) X 10~^ cm~^ for the para line at 269.3 ;um while the ortho 
line at 538.3 /xm has Wo = (0 it 6) x 10~^ cm~^. A given equivalent width in 
the ortho line implies 4 times more column density than the same equivalent 
width in the para line because the frequency is half as large and the transition 
matrix element l/ijjp is half as large. Thus ortho:para column density ratios 
ranging from to more than 3:1 are all consistent with these equivalent widths. 
If both lines are optically thick, as seems likely, then Wo should be one half of Wp 
because of the frequency factor, which is also consistent with the observations. 
In the optically thick case we have no information on the ortho:para ratio. 

Dr. Dwek: How do you interpret the two-temperature dust fit in terms of a 
physical dust model? 

Prof. Wright: It could be caused by a shoulder in the dust emissivity law. This 
would explain the constancy of the cold dust temperature and the cold/warm 
opacity ratio. 
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